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ABSTRACT 

Aiming to perforin a study of the warm dust and gas in the luminous blue variable star 
G79.29+0.46 and its associated nebula, we present infrared Spitzer imaging and spectroscopy, 
and new CO J = 2 — )■ 1 and 4 — >■ 3 maps obtained with the IRAM 30m radio telescope and with 
the Submillimeter Telescope, respectively. We have analyzed the nebula detecting multiple shells 
of dust and gas connected to the star. Using Infrared Spectrograph-Spitzer spectra, we have 
compared the properties of the central object, the nebula, and their surroundings. These spectra 
show a rich variety of solid-state features (amorphous silicates, polycyclic aromatic hydrocarbons, 
and CO2 ices) and narrow emission lines, superimposed on a thermal continuum. We have also 
analyzed the physical conditions of the nebula, which point to the existence of a photo-dissociation 
region. 



Subject headings: 
lines and bands 

Introduction 



circumstellar matter - dust, extinction - ISM: individual objects (G79.29+0.46) - ISM: 
- stars: evolution - stars: mass-loss 



Luminous blue variable (LBV) stars are mas- 
sive objects which, as they evolve from the main 
sequence, undergo a short period of extremely 
high mass loss (up to 10"'^ yr~^), sometimes 
accompanied by so-called giant eruptions, such 
as the well-known nine teenth-century outburst of 
77 Car (jSmith fc Owoc ki 2006). This mass-loss 
strongly influences the further stellar evolution 
and leads to the formation of extended circumstel- 
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lar nebul a (see Humphrevs fc Davidsonl[l994 , for a 
review). IClark et al.l (|2005h reported a total of 12 
galactic LBV stars and 23 LBV candidates. From 
the whole sample, however, only 20 have nebular 
circumstellar emission. Such nebulae have both 
gaseous and dusty components, and two kinds of 
geometry, spherical and bipolar, have been recog- 
nized. 

Since LBV nebulae are strong emitters in the 
mid- and far-infrared, the physical and chemical 
conditions of the material surrounding LBV stars 
can be traced studying their spectral characteris- 
tics in this wavelength range. Thus, a few of these 
nebulae were spectroscopically studied with the 
Infrared Space Observatory {ISO) satellite, reveal- 
ing the presence of solid-state features character- 
istic of polycyclic aromatic hydrocarbons (PAHs) 
and silicate dus t in both forms, arnqrphous and 



crystalline ( e.g. iLamers et al.lll996a[ Voors et aD 
I999I 2000bf) . Together with these, numerous for- 
bidden lines on top of a sm ooth continuum were 
detected (e.g. Lamers et al.l[l99 6b'). 

With the adven t of Spitzer Space Telescope 
( Werner et al. 20041 ) , unprecedented high sensitiv- 
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ity and resolution imaging and spectroscopy have 
been possible in the infrared, making feasible the 
study of LBV nebulae with much great detail. 
Thus, recently new works have been published 



(iMorris. P .. et al."2008': 'Smithl|2007: lUmana et al 



2009t iGva raniadzc ct al. 2009) based on these new 
data, which have represented an important step 
forward in the knowledge of this brief evolution- 
ary phase of high-mass stars. 

The molecular gas in LBV nebulae is also a 
subject of interest. It may trace radiatively af- 
fected and/or shocked regions, which may tell us 
about the evolution of the progenitor star. In 
addition, the chemical evolution and the forma- 
tion of circumstellar molecular gas may also be 
determined by molecular studies. Gaseous CO 
was firstly detected around the AG Car nebula 
(INota et al. 20021). and aft er then in some other 
objects (jRizzo et al.ll2008a[l . Highly excited am- 
monia was also detected in the Homunculus Neb- 



ula (,Sniith et al.ll2006i) 



Concerning the molecular count erpart, the 
most studied source is G79. 29-1-0. 46 (Rizzo et al 



2008bl hereafter Paper I). In this object, the de- 
tection of the CO J = 3— !>2 line surrounding the 
nebula unveiled the presence of moderately dense 
(lO-^-lO^cm-^) CO gas, probably affected by a 
C-shock at a velocity of ^^ISkms^^. This CO- 
ring-like structure has a diameter of ^5'.2, and it 
is strikingly homogeneous in shape. 

The origin of the nebula of G79.29-f 0.46 is not 
clear. From Very Large Array ( VLA) obser v ations 



of radio continuum and H I, iHiggs et al.l (|l994l ) 



showed the thermal nature of the continuum emis- 
sion and suggested that the nebula is an ionized 
shell of ~ 15 Mq of swept-up interstellar mate- 
rial. However, after examining th e IRAS high- 



resolu tion images of G79.29+0.46, IWaters et al 



(Il996l) concluded that it is a detached shell due to 



an epoch of high mass loss ('^5xl0~ M0yr^ ) 
during the red supergiant or LBV phase, followed 
by a less intense mass-loss period. 

Concerning the central object, it has been 
spectroscopically studied in the o ptical and in 



the near-infrared wavelen gth range (jWaters et al 



1996t IVoors et al.ll2000b[ ). Its LBV classification 
is supported by the high luminosity (>1O^L0), 
the moderate stellar wind (94 -110 km s~^), and 
the current high m ass-loss rate (10~^ yr"-*^) 
( Waters et al. 19961 ). Unfortunately, its effective 



temperature {T^ff) is not well-known because it 
is highly model dependent. The detection of a 



He I emission line near Bra ( Waters et al. 19961 ) 



suggests a Teff well above 10,000 K. On the other 
hand, the LBV nature of the central star may 
limit the Teff to a maximum of 30,000 K. 

In Paper I, a careful discussion about the dis- 
tance to G79. 294-0. 46 was done. We will assume 
a distance of 1.7kpc throughout all the paper. 

The fingerprint of stellar evolution onto the 
circumstellar/interstellar medium (ISM) can pro- 
vide a valuable input to learn about the processes, 
which govern the star itself and also affect the evo- 
lution of the surrounding gas and dust. The aim 
of this paper is to trace the distribution of warm 
dust, PAHs, and gaseous CO in relation with 
G79. 29-1-0. 46 and its associated nebula. To do 
this, we present infrared Spitzer imaging and spec- 
troscopy, and new submillimeter CO J = 2 — > 1 
and 4 — )• 3 maps obtained with the IRAM 30m 
radio telescope and the Submillimeter Telescope 
(SMT), respectively. In comparison to Paper I, 
this new data set improves the angular resolution 
and the analysis of the physical conditions. 

2. Spitzer data 

The Spitzer data used in this paper have been 
obtained from the public Spitzer archiv^. 

2.1. Imaging 

Two cameras are on board of Spitzer: the In- 
fraRed Array Camera (IRAC; iFazio et al.l \2004 ) 
and the Multiba nd Imaging Photometer (MIPS; 
Rieke et al. which provide imaging capabil- 

ity from 3 to 200 /im. 

Near-infrared images were obtained on 2007 
July 05 using the four IRAC bands at 3.6, 4.5, 
5.8, and 8.0 ^m, as part of the program 30188 
(AOR 17330688). Mid-infrared images at 24 and 
70 ^m were also obtained using MIPS as part of 
the program 40184 (AOR 22510592) on 2007 De- 
cember 02. We have worked with the standard 
post-BCD data sets, for both IRAC and MIPS im- 
ages, as retrieved with Leopard from the Spitzer 
archive. They were processed with the version 
S16.1.0 and S16.1.1, respectively. Aladird was 



^jhttp: //archive ■ spitzer ■ caltech. edu] 
• ^http : //aladin . u- strasbg ■ f r] 
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used to cut the images to show a field of view 
of 9' X 9' centered on the LBV star. These images 
are shown in Figure [TJ The brightness scale of ev- 
ery filter has been individually selected to better 
investigate the nebular emission of G79.29+0.46. 

2.2. IRS Spectra 

The Spitzer satellite has a third inst rument on 
board , the Infrared Spectrograph fIRS: Houck et al 
20041 ) ■ It provides crucial spectroscopic informa- 
tion in the wavelength range from 5 to 40 /xm. 

Three targets in the field of G79. 29-^0.46 have 
been spectroscopically observed in four observa- 
tional programs: a target including the exciting 
star and its surroundings (hereafter called 'the 
central object'), a target centered at a position 
where the ring nebula is most prominent (here- 
after called 'the shell'), and a nearby infrared dark 
cloud (hereafter called 'the IRDC). These obser- 
vations are tabulated in Table[Tl and the positions 
of the IRS slits are shown in Figure [2j It should be 
noted that all three, the central object, the shell, 
and the IRDC, include emission from nearby re- 
gions than the targets themselves. 

The spectra were obtained by using all pos- 
sible modules: short-low (SL), short-high (SH), 
long-low (LL), and long-high (LH). There are low- 
resolution spectra (SL and LL) for the three stud- 
ied targets, which cover the spectral range from 
5.2 to 38.0 /xm. However, high-resolution spectra 
(SH and LH), covering the spectral range from 
9.2 to 37.2 /im, are only available for the central 
object. We have worked with the standard post- 
BCD data sets. We have cleaned the resulting one- 
dimensional spectra for residual bad pixels flagged 
by the pipeline and visually inspected to identify 
spurious jumps and glitches. No spurious jumps 
were detected and the glitches found were interpo- 
lated. Finally, we have merged the different orders 
of each module into one final spectrum per mod- 
ule using an IDL code developed to this end. We 
have not corrected for basic aperture and slit loss. 
Spectra are shown in Figures |3]-[5l 

In the case of the central object, the target 
has been observed on three different occasions 
(see TablelJ). The first time (AOR 7557888), the 
IRS slits were not correctly positioned over the 
central object. The second time (AOR 9761024), 
there should have been a flux calibration prob- 



lem, since the two orders do not match. The 
third time (AOR 17333248), the data do not show 
any problem. Thus, we have used these last spec- 
tra for the analysis, corresponding with the mod- 
ules SL, SH, and LH. However, on this last oc- 
casion LL spectrum was not taken, so we have 
used for the analysis the LL spectrum correspond- 
ing to the AOR 9761024. The SL spectrum of the 
AOR 9761024 and those of the AOR 7557888 have 
been used to conflrm the detection of the solid 
spectral features and the emission lines (see Sec- 
tion HJl). 



3. CO observations 

3.1. CO J — 2^1 Line Emission 

We used the IRAM 30m radio telescope at Pico 
Veleta (Spain) in 2008 February, to map the CO 
J = 2 ^ 1 line emission in a region of ^ 7' x 7' 
around G79. 29-1-0. 46. The front end wa.s the nine- 
beam HERA receiver ( Schuster et al.l |2004|) . In 



order to optimize both observing time and sam- 
pling, we used the on-the-fly mode. A total of 
four individual maps, in orthogonal sampling di- 
rections, were obtained and later combined in a 
single data cube using the GILDAS package. The 
VESPA autocorrelator was used as back end, in 
frequency switching mode. The flnal bandwidth 
was 44 MHz (equivalent to ^ 57kms~^ at the rest 
frequency) , and the spectral resolution was 80 kHz 
(^O.lkms^^). Due to the frequency switching, 
mesospheric CO emission had to be removed from 
all the spectra. After folding, low-order polynomi- 
als were fltted to remove baselines using CLASS 
software package (also included in GILDAS). Typ- 
ical rms was 0.2 K, in scale. 

Calibration was performed during the obser- 
vations using two absorbers at different temper- 
atures and later checked u p on standard sources 
( Mauersberger et al. 19891 ). The uncertainty in 
calibration is below 15%. The angular resolution 
of the radio telescope at 230 GHz half-power beam 
width (HPBW) is 11". However, we have further 
degraded the final map down to 18", in order to 
reduce the map noise and to get the same HPBW 
corresponding to the CO J = 4 3 map. 

Some other emission lines were observed during 
the same observing session. The results from these 
lines will be released in a future paper (J. R. Rizzo 
et al. 2010, in preparation), together with a more 
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Fig. 1. — Spitzer images of G79.29+0.46. Wavelengths are indicated at the upper part of each image, and 
angular offsets are referred to the central star, at (R.A., decl.)j2ooo = (20'*31'"42!3, +40°21'58"). A linear 

brightness scale has been artificially selected for each individual image to better show weak and diffuse 
features, instead of point-like objects. Note an outer thin second shell clearly visible only at 24 /xm. 




Fig. 2. — Overplot on the MIPS 24 image of the IRS spectral shts for all modules, as shown by Leopard 
software. The field of view is 14'xl4' centered on the LBV star. Only the observations used in the analysis 
are shown. The three different positions correspond with the central object, the shell, and the IRDC. 



5 



Table 1 

irs observations in the field of g79.29+0.46 



ID 


AOR 


Date 


R.A 




Decl 






Target 


Module 


A Range 


Slit Size 


Pipeline 


program 






J(20Q0) 


,1(2000) 








(^im) 


(aresee) 




666 


7557888 


2003 Nov 14 


20*^ 31™ 


42=0 


+40° 21' 


42' 


'0 


The central object 


SL 


5.2-14.5 


3.6x57 


S15.3.0 




















SH 


9.9-19.6 


4.7x11.3 


S15.3.0 




















LL 


14.0-38.0 


10.5x168 


S17.2.0 


1407 


9761024 


2004 May 11 




42! 1 


+40° 21' 


58' 


'8 


The central object 


SL 


5.2-14.5 


3.6x57 


S15.3.0 




















LL 


14.0-38.0 


10.5x168 


S17.2.0 




9761280 


2004 May 14 




48! 9 


+40° 23' 


15' 


'4 


The shell 


SL 


5.2-14.5 


3.6x57 


S15.3.0 




















LL 


14.0-38.0 


10.5x168 


S17.2.0 


3121 


12096768 


2004 Dec 07 


20** 31™ 


45! 6 


+40° 18' 


43^ 


'6 


The IRDC 


SL 


5.2-14.5 


3.6x57 


S15.3.0 




















LL 


14.0-38.0 


10.5x168 


S17.2.0 


30188 


17333248 


2006 Jun 28 


20*^ 31™ 


42! 5 


+40° 21' 


56' 


'0 


The central object 


SL 


5.2-14.5 


3.6x57 


S15.3.0 




















SH 


9.9-19.6 


4.7x11.3 


S15.3.0 




















LH 


18.7-37.2 


11.1x22.3 


S17.2.0 



complete analysis of the CO J = 2 — ;> 1 results. 

3.2. CO J — 4^3 Line Emission 

The observations of the CO J = 4 — > 3 line were 
made in 2006 June with the SMT, located on 
Mount Graham, Arizona. The telescope is a 
10.4 m diameter primary with a nutating sec- 
ondary. Design, optics, and structural aspects 
of the teles c ope w ere extensively described by 



iBaars et al 



(Il999l) . The single- channel receiver 
SORAL-490 was used under excellent weather 
conditions {t4qichz '-^0.4-0.5). As the back end, 
we used one of the acousto-optical spectrometers 
(AOSCs), which provides an usable bandwidth of 
about 250 MHz (corresponding to 162 km s^^ of 
coverage), and a spectral resolution of 367kHz 
(0.24kms~i). 

The receiver was calibrated with ambient tem- 
perature and a liquid nitrogen-cooled absorbing 
load. Calibration of the final data has uncertain- 
ties below 20%. System temperature was mea- 
sured with the st andard sky/amb i ent t empera- 
ture load method ("K utner fc UhchI Il98ll ). Cali- 
brations were performed every 15-20 minutes to 
track variations in atmospheric opacity. Typical 
system temperature during the observations was 
in the range of 1000- 1200 K, depending on the 
elevation and the atmospheric opacity. 

The pointing was checked within 2 hr on a 
planet or other bright continuum source. Correc- 
tions to the pointing model were not more than 
4" in either azimuth or elevation. Planet obser- 
vations were also used to compute the main-beam 



efficiency at the observed frequency. 

Only a small area toward the southwest of 
G79. 29-1-0. 46 was observed, where the most in- 
teresting CO features appear (Paper I). Position 
switching was employed, using the same reference 
position for all observed points. Due to time re- 
striction, the CO J — 4^3 map was sampled only 
by one beam (18") in both equatorial coordinates, 
and hence it is not fully sampled. The reduction 
processes were similar to those of IRAM's data, 
already described in Section [3Tl 

4. Results 

4.1. Infrared images 

The distance assumed for G79. 29-1-0. 46 is con- 
sistent with its probable identification as a mem- 
ber of Cyg 0B2. Located at the southern part 
and close to the border of Cyg 0B2, G79. 29-^0.46 
is behind the Great Cygnus Rift, just to the 
northwest of the intense radio source DR 15, 
which causes high visual extinction from 4 to 
at least 10 mag. Thus, the infrared emission of 
G79. 29-1-0. 46 is expected to be heavily contami- 
nated by foreground/background emission, as it is 
clearly noticeable in all Spitzer images. In addi- 
tion, a large-scale infrared emission is also present 
toward the east of the images. This large-scale 
emission increases its relative intensity with re- 
spect to the background with the wavelength, and 
is similar to the large-scale CO gas emission noted 
in Paper I at millimeter and submillimeter wave- 
lengths. It was r elated to the cloud around DR 15 
(jOka et al.ll2001l ). Due to this strong contamina- 
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Fig. 3. — IRS low-resolution spectra (modules SL 
and LL) of the three studied targets. The most 
prominent spectral lines and solid features are la- 
beled in the figure. Upper panel: the spectra of 
the central object after subtracting the LL spec- 
trum of the shell (see the text). The dashed line 
corresponds to the original LL spectrum. The 
dotted (blue) lines correspond to the emission of 
two blackbodies, one at 20,000 K and the other at 
925 K. The dash-dotted (red) line is the combina- 
tion of these two blackbodies. Middle panel: the 
spectra of the shell. The dotted (blue) lines corre- 
spond to the emission of two blackbodies, one at 
65 K and the other at 132 K. The dash-dotted (red) 
line is the combination of these two blackbodies. 
Lower panel: the spectra of the IRDC. The dot- 
ted (blue) lines correspond to the emission of two 
blackbodies, one at 66 K and the other at 187 K. 
The dash-dotted (red) line is the combination of 
these two blackbodies. 
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Fig. 5. — IRS high-resolution spectra (modules SH 
and LH) of the central object. The most promi- 
nent spectral lines are labeled in the figure. 
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tion, we refrain from computing the dust mass, 
as it is a hard task and the results might not be 
rehable enough. 

The central object is detected in all Spitzer 
bands but MIPS 70 /im. We used MOPEX sup- 
ported software to perform point-spread function 
(PSF) photometry in the MIPS 24 /im image, ob- 
taining a flux density of 0.89 Jy for the central ob- 
ject. Due to saturation problems, it was not pos- 
sible to obtain the fluxes of the central object at 
shorter wavelengths. 

The ring nebula previously detected by the 
IRAS an d Midcourse Space Experiment (M SX) 
satellites (|Waters et al.lll996t lEgan et al.|[2002D is 
discernible from 3.6 to 70 ^m. The nebular emis- 
sion dominates the image at 24 /im, being just 
above the background level at shorter and longer 
wavelengths. The nebula presents a pretty spheri- 
cal shape with a diameter of ~ 3.'7 (~ 1.8 pc). The 
shell looks like incomplete to the south, probably 
due to a foreground IRDC, which is known to pos- 
sess a very large visual extinction ( Redman et al.l 
20031 ). However, it is not clear whether this IRDC 
is part of the large complex incl uding DR 15 or it i s 
an unrelated foreground cloud ( Egan et al. 1998[) . 

Thanks to the high angular resolution and sen- 
sitivity of these Spitzer images, we have discovered 
an outstanding new, outer shell with almost spher- 
ical shape as well. This outer shell is only visible 
in the 24 /im image (Figure [T]) and has a diameter 
of ~6f8, (^3.4 pc). During the proces s of writ- 
ing this paper. iGvaramadze et al. ( 20091) reported 
a diffuse halo associated with G79. 29-1-0. 46, which 
corresponds to the inside of this newly discovered 
shell. 

4.2. CO maps 

As we show below, the CO data shown here 
improve the results gathered in Paper I, due to 
a higher angular resolution and lower rms noise 
compared to those in Paper I, providing a higher 
signal-to-noise ratio (S/N). In addition, mapping 
the CO J = 4 3 line emission allows us to test 
a wider range of excitation conditions, due to 
the higher energy levels involved and the criti- 
cal density needed to excite this line. Accord- 
ing to the aim of this work, here we provide the 
most relevant morphological features, postponing 
a more detailed dynamical study to a future paper 



(J. R. Rizzo et al. 2010, in preparation). 

The whole velocity range of emission in the 
J = 2— T'l line goes from -25 to 15kms~^. The 
most intense emission corresponds to the Local 
Arm, where the confusion effects are presumably 
the highest. The most outstanding features in re- 
lation to G79.29+0.46 appear in velocity ranges 
similar to those shown in Paper I (see Figure 4 
of that paper). Figures IH] and [7] show the CO 
J = 2 — > 1 and 4 — > 3 distributions in the veloc- 
ity ranges from -3.2 to -2.6 and from -1-3.4 to 
-1-4.0 kms~\ respectively. The CO J = 2^1 line 
contours are superposed to an image at 24 /im in 
gray scale. 

For the J = 2 — > 1 data, a Gaussian convolu- 
tion has been applied, down to a resolution of 
18". This convolution reduces the noise map and 
equals the angular resolution of the 7 = 4^-3 line. 
In both figures, the inset in the J = 2 — ^ 1 maps 
corresponds exactly to the mapped area in the 
J = 4 — )• 3 line, shown on the right. 

In Figure [6l an arc-like feature dominates the 
map, in good coincidence with the 24/Ltm inner 
shell, especially in the western half. The J = 4 — 3 
line emission depicts part of a remarkable thin 
shell, unresolved even at 18" of resolution. 

The shocked clump, unresolved in Paper I, is 
clearly seen in Figure [T] In this velocity range, 
all the emission is outside the 24 /im inner shell. 
The southwestern feature, also observed in the 
J = 4 — >■ 3 line emission, looks like part of a sec- 
ond CO shell. Note that this feature is located 
outside the shell shown in Figure |6l 

The emission at velocities intermediate to those 
of Figures [S] and [7] does not show any shell-like 
pattern, and consequently we think that the two 
shells are different features, and not the conse- 
quence of a mere velocity gradient. Summarizing, 
at the southwestern part of the field, we see two 
CO concentric slabs, which fill the space between 
the two shells detected at 24 /im. 

4.3. Infrared spectra 

The infrared spectral data show a rich variety 
of solid-state features and narrow emission lines, 
superimposed on a thermal continuum. 

Using the low-resolution IRS spectra, it is pos- 
sible to study the spectral energy distribution 
(SED) of the three observed targets (Figure [3]). 
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C0(3-l) 




Fig. 6.— Left: CO J = 2^1 map in the surroundings of G79.29+0.46, as obtained with the IRAM 30m 
radio telescope. The map is centered on the LBV star, and the velocity range of integration is from -3.2 
to -2.6 km s"'^. The pointed square indicates the observed field, and gray scale corresponds to the 24 /im 
emission. Contours are 40%-90% of the peak value (11. 8K in scale), in steps of 10%. Right: CO 
J = 4— >-3 map in the same velocity range, corresponding exactly to the inset in the left panel, as obtained 
with the SMT 10 m radio telescope. The inset has a size of 80" x 116". Contours are 25 % -95 % of the peak 
value (1.8 K in scale), in steps of 10%. In both maps, the angular resolutions are shown as white circles 
in the lower right. It is striking the high correlation of the CO emission in this velocity range with the 
infrared inner shell. The CO shell is especially thin in the J = 4 — >^ 3 map, unresolved by the beam. 
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Fig. 7. — Same as Figure [6l but in the velocity range from +3.4 to +4.0 km s^^. Left: contours are 30%- 
90% of the peak value (10.8 K in scale), in steps of 10%. Right: contours are 25% -95% of the peak 
value (3.4 K in scale), in steps of 10%. This velocity range corresponds to a shocked clump reported in 
Paper I. After increasing the angular resolution, a second, outer shell is clearly visible in both lines. 
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The SED of the central object (upper panel) shows 
that the flux decreases from 5 to about 19 fj,m and 
increases again to longer wavelengths. This be- 
havior may be roughly explained by two emitters 
at two different temperatures, i.e., hotter (shorter 
wavelengths) and colder material (longer wave- 
lengths). The middle and lower panels of Fig- 
ure [3] correspond to the spectra of the shell and 
the IRDC, respectively. These two last SEDs are 
dominated by the colder component. 

The three SEDs present an astonishing similar- 
ity for A > 25 /im. This suggests that this emission 
is due to interstellar material which is present in 
all regions and extends beyond the observed area. 

Figure [2] shows the position of the IRS spectra 
slits for the three regions observed. It is clear that 
the SL and LL spectra of the central object are 
affected by the shell emission. The effect is maxi- 
mum at longer wavelengths, since the emission of 
the shell at shorter wavelengths is negligible when 
compared with the emission of the central object 
(see Figure [3]). Thus, in order to remove the con- 
tribution of the shell emission to the central ob- 
ject SED, we have scaled the LL spectrum of the 
shell to fit that of the central object at the longest 
wavelengths (>32/im) and later subtracted from 
the last one. The result is shown in Figure [3] (up- 
per panel) with a solid line. The subtracted cen- 
tral object SED shows the decreasing tail of the 
central object emission. In addition, several spec- 
tral lines ([Ne III] at 15.5 fim, [S III] at 18.71 and 
33.4 /xm, H2 at 28.2 ^m, and [S III] at 34.8 ^m) 
have disappeared of the central object spectrum. 

4-3.1. Solid-state Features 

The SEDs of the central object and the IRDC 
reveal an absorption feature at around 9.7 /xm, 
more prominent in the case of the central object, 
particularly clear after subtracting the shell con- 
tribution (Figures [3] and S]) . It corresponds to the 
well-known 9.7 /xm amorphous silicates feature due 
to the Si-0 stretching in the dust grains. The 
presence of amorphous silicates was also reported 
in th e dusty nebula surrounding the LBV star s HR 
Car (jLamers et al."l996a': ' Umana et aLll2009f ) and 
Wra 751 (Voors ct al . 2Q003) . Note , however , that 
these stars present the feature in emission instead 
of absorption as in the case of G79.29-f 0.46, which 
is puzzling and has to be explained. In addition, 
the SED of the shell seems to present a shallow 



depression at around 9.7 /im that could be inter- 
preted as a weak absorption due to amorphous 
silicate dust. 

The broad, strong emission features at 6.2, 7.7, 
8.6, 11.3, and 12.7/xm present in the SL spec- 
trum of the shell (See Figure lU upper panel) 
are due to vibrational modes of free-flying PAH s 
(|Leger fc Pugetl 11984 lAllamandola et all [l985l) . 
The same features seem to appear in the spec- 
trum of the IRDC, but the flux is too low for 
reliable detections. These PAHs features suggest 
the presence of very small C rich dust grains or 
molecules in the shell of G79. 29+0.46. PAHs have 
also been dete cted in the shell of the galactic LBV 
stars AG Car (IVoors et all bOOOal) and HR 168625 
(|Skinneij[l997l ). and in t he LBV star R71 of the 
Large Magellanic Cloud (jVoors et al. I I1999I) . 

Another interesting solid-state feature, present 
only in the spectrum of the IRDC, is the CO2 
bending mode at about 15 /im (see Figure |31 lower 
panel). Since the formation and evolution of in- 
terstellar ices is strongly dependent on local con- 
ditions, the presence of CO2 ice in the spectrum 
of the IRDC indicates that the physical condi- 
tions are rather different from those in the cen- 
tral object and in the shell. CO2 ice has never 
been detected in the surroundings of a massive 
evolved star, but is usual in other astronomical ob- 
jects, such as high- and intermediate-mass proto- 



jects, sucn as nign- ana intermediate-mass proto- 
stars (lEhrenfreund et alll998HPeeters et al.112002 : 



Lefloch et al. 20081 ). The CO2 absorption feature 
found could be coming either from the quiescent 
regions of the IRDC (south of the LBV nebula) or 
from the compact source lying within the south- 
ern slit (see Figure [2]). Note that this compact 
source was reported as a young stellar object by 



Redman et al.l 



fortea a 
(I2OO3I) . 



4.3.2. Emission lines 

The IRS high-resolution spectra of the central 
object are full of fine structure lines that may pro- 
vide important constraints to the physical condi- 
tions in the region. The most prominent ones are 
also detected in the low-resolution spectra and are 
labeled in the figures. For each line, a first-order 
baseline has been defined and the line has been 
fitted to a Gaussian. To do thi s, we have u sed an 
iterative software developed bv lBavol(|2009l) . The 
identification and the measured fluxes are summa- 
rized in Table[21 including possible multiple iden- 
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tifications. Flux errors are between brackets. All 
line fluxes for the central object have been mea- 
sured in the high-resolution IRS spectra but H I 
at 7.46 /im. In the case of the shell and the IRDC, 
line fluxes have been measured in the respective 
low-resolution IRS spectra. There is a mismatch 
between the SH and LH spectra of the central ob- 
ject which is usually attributed to a difference in 
the two module apertures, and which has been 
taken into account when comparing fluxes of SH 
with those of LH spectra. 

The collection of spectral lines found is sim- 
ilar to that shown by the LBV star HRCar 
(jUmana et al."20 09|) and the LBV candidate HDE 
316285 ( |Morris. R. et all |2008|) . The nebular 
spectra of HDE 316285 show, for the first time 
in this kind of objects, pure rotational lines of H2 
arising from optically thin quadrupole transitions, 
from S(0) to S(6). We have also clearly detected 
the two lower rotational lines, S(0) and S(l), in 
the spectra of the central object and the shell, and 
S(0) in the spectrum of the IRDC. Although these 
lines may arise from a foreground component, or 
even DR15, the presence of H2 in the shell or in 
the circumstellar material is a possibility which 
has to be explored in follow-up observations. 

5. Discussion 

5.1. Multiple shells in G79.29+0.46 

The data shown in this paper clearly show mul- 
tiple dust and gas shells connected to the cen- 
tral star of G79.29-h0.46. The MIPS 24 /im image 
shows two well-defined dust shells. Figure [6] shows 
how the CO shell coexists next to the inner dusty 
shell, while Figure [7] shows another slab located 
between the two dust shells. 

To explore the possibility of the existence of a 
more inner shell closer to the LBV star, we have 
roughly fitted the SED of the central object af- 
ter subtracting the contribution of the shell with 
a combination of two blackbodies. Two ranges of 
temperatures were tried, from 10,000 to 30,000 K 
for the stellar contribution and from 200 to 1500 K 
for a possible inner envelope. Fitting does not vary 
significantly in the considered ranges of tempera- 
ture corresponding to the hotter blackbody. How- 
ever, the best fit regardless the stellar temperature 
was always obtained with a second blackbody of 
~900K. This result is compatible with the exis- 



tence of an inner shell unresolved by the Spitzer 
images. Figure |3] shows the best fit for a stellar 
temperature of 20,000 K. This conclusion is sup- 
ported by the detection of H2 lines in the spectra 
of the central object. It is necessary a thick dust 
envelope to be able to explain the existence of H2 
so close to the LBV star. This thick dust envelope 
would shield the molecule gas from the strong UV 
field of the star. 

From the morphology, it seems that the two 
dusty shells resolved by the Spitzer images are in- 
teracting with the surrounding medium (see Fig- 
ure [1]), presumably with material belonging to the 
DR15 cloud. Figure [7] also shows clear evidence 
of the erosion of the ISM by the stellar wind, spe- 
cially at the southwestern region. At this region, 
we already reported in Paper I the presence of a 
shock front. 

In Paper I, we estimated an age of the nebula 
between 10'^ and 10^ yr. This is in well agreement 
with the high sphericity shown by the two dust 
shells. Thus, the mass-loss events which have pro- 
duced the multiple shells have occurred relatively 
recently, i.e., in the current evolutionary stage. 

In order to obtain a dust temperature, we have 
also fitted the continuum emission of the shell and 
the IRDC, with a combination of two blackbod- 
ies at different temperatures. In the case of the 
shell, a blackbody of 65 K and another of 132 K 
were necessary. For the IRDC, one of 66 K and 
another of 187 K were used. The resulting fits 
are shown in Figure [3l In both cases, two dust 
components at different temperatures were found, 
one common component at 65 K and another 
hotter. Probably, the cold component is due to 
foreground interstellar material belonging to the 
H II region DR15, as typical dust temperatures 



of H II regions are around 60 - 70 K (jPeeters et al 



20021 ). This is the main dust component and would 



explain why both spectra have the same overall 
shape at longer wavelengths. Regarding the hot 
component of the shell, its temperature is simi- 
lar (same order) to the temperatures e stimated for 



the d ust in other LBV nebulae (300 K: IVoors et al 



I999I ) and in supernova rem nants (SNRs) asso- 
ciated with shoc ks (^ llOK; iTappe et all l2006t 
2OO9I) . Note that the temper- 



Ghavamian et al 



ature of the hot component of the shell is slightly 
lower than the temperature of the hot component 
of the IRDC. Although these temperatures should 
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Table 2 

Line Identification and Measured Fluxes 



Line 


y 




Flux (10-14 erg 






(//m) 


The ( 'ciilral Object 


Tlic- SIk-U 


TIk- liU)C 


r A TTl Fm TTl 

[Ar iij; [Ne ii] 


6.97 




7.7 (0.6) 




H I 6-5 8-6 


7.46 


960 (90) 






[Ar III] 


8.99 


6 (2) 




He II 10-9 


9.67 






2.6 (0.8) 


He II 24-16; H I 12-8; [S IV] 


10.50 


19 (3) 






[Ni II] 


10.68 


7.7 (0.6) 






riVT TTTl 

[Ne 111] 


10.88 


16 (2) 






He II 23-16 18-14; [Ni I]; H I 9-7 


11.30 


28 (2) 


58 (8)'' 




H I 15-9 


11.53 


8 (2) 






He II 14-12; H I 7-6 


12.37 


77 (7) 




1.9 (0.2) 


[Ne II] 


12.81 


18 (2) 


98 (8)'' 


4.6 (0.2)'' 


H I 17-10; H I 11-24 


13.93 


6.7 (1.0) 






H I 13-9 


14.17 


5.1 (0.6) 






H I 23-11 


14.29 


3.0 (1.0) 






H I 16-10 


14.94 


6.7 (1.3) 






[Nc III] 


15.55 


13 (2) 


29 (3) 




He II 20-16- H I 10-8 


16.20 


10 (2) 






He II 24-18; H I 12-9 


16.87 


3.6 (0.3) 






H2 0-0 S(l) 


17.03 


4.0 (0.5) 


19 (6) 




H I 18-11 


17.60 


6.7 (2) 






[Fc II] 


17.90 


8.4 (0.4) 






SIII 


18.71 


24.0 (1.0) 


184 (15) 


18 (4) 


He II 23-18 16-14; H I 8-7 


19.05 


24 (5) 






H I 20-12 


20.52 


4.5 (0.6) 






[Ar III] 


21.82 


5.2 (0.6) 






He II 17-15; H I 13-10 11-9 


22.32 


27 (2) 






H I 15-11 


23.86 


27.0 (1.4) 






[FI] 


24.75 


9.0 (0.5) 






H I 21-13 


24.98 


4.1 (0.8) 






[Fell] 


25.97 


23 (3) 






H I 24-14 


27.08 


7.3 (0.5) 






He II 18-16; H I 9-8 


27.79 


29 (2) 






Ha 0-0 S(0) 


28.21 


13 (4) 


25 (6) 


11 (2) 


H I 14-11 


28.85 


7.7 (0.8) 






H I 19-13 


28.97 


14 (3) 






He II 24-20; H I 12-10 


29.83 


10.3 (0.7) 






[S III] 


33.45 


78 (11) 


740 (40) 


229 (15) 


[Sill] 


34.80 


41.5 (1.1) 


150 (40) 


66 (6) 



^Pcak wavelength measured on the stellar spectrum. 
''Flux measured affected by PAHs emission. 
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be seen with caution due to the presence of the 
CO2 ice absorption feature, the shghtly higher 
temperature in the IRDC could be due to the pres- 
ence of a young stellar object in the slit (see Fig- 
ure H]). 

The most external dust shell is only visible in 
the 24 fim image. It is expected that the dust in 
this shell would have a lower temperature than the 
dust in the inner shells. Thus, the dust tempera- 
ture of the most external shell should not be well 
above 120 K. Assuming that the mid-infrared im- 
ages are tracing the thermal continuum emission of 
the dust, this would explain why it is not detected 
at shorter wavelengths. The fact that this shell is 
neither detected at 70 could be due to a combi- 
nation of strong foreground/background emission 
at this wavelength, lower sensitivity of MIPS in 
this band, and a non-appropriate data reduction 
of the image. 

The 9.7 /^m amorphous silicates feature in ab- 
sorption is clearly present in the SL spectra of 
the central object and the IRDC, and is tenta- 
tively detected in the shell. Amorphous silicates 
in LBV stars are supposed to be formed at tem- 
peratures well below the silicates melting point in 
the outer layers of the star which cool very rapidly 
while the material is being swept up by a fast 
wind. The silicates would indicate that the stellar 
ejecta was oxygen rich (C/0 < 1), as expected for 
LBV stars which are sup posed to have eje cta en- 
riched by the CNO cycle (|Nota et al]|l997t ). How- 
ever, G79. 29-1-0.46 shows the 9.7 /im amorphous 
silicates feature in absorption instead of in emis- 
sion as the others LBV stars do. A plausible ex- 
planation is that the absorption of amorphous sil- 
icates is due to the foreground ISM, as was re- 
cently found in the c ase of the SNR N157B by 



Micelotta et al 



(|2009l ). This would explain the 
presence of this feature in absorption in the IRDC 
as well. However, another possible explanation for 
the absorption shown in the central object spec- 
trum is the presence of amorphous silicate dust in 
an inner nebula unresolved by the Spitzer images. 
The presence of crystalline silicates in LBV neb- 



ulae seems to be a general trend (jWaters et al 
19971 ). This is suggestive of the grains being 



formed in hi gh-density and lo w-temperature en 



vironments ( Voors et al. I1999I) . Therefore, and 



silicates is explained if the dust is formed in high- 
density, slow-expanding environments close to the 
star, rather than in the swept-up shell. This kind 
of env ironment is typical of the red supergiant 
phase ([Waters et al.lll998( ). for which wind veloc- 
ities are between 20 and 30kms~^. In the case 
of G79.29-h0.46, and similarly to the LBV can- 



didate HDE 316235 (jMorris. P.. et al 2008) and 



the L BV stars Carina and HR Car (jUmana et al 



20091) . its spectrum does not show any signature 
of crystalline silicates. The lack of crystalline sil- 
icates could be explained if the nebular material 
was rapidly ejected and could not crystallize un- 
der conditions of rapid c ooling and low monomer 
densities, as suggested bv lMorris. P.. et"aL I (l2008h 
for HDE 316235. Another possible explanation 
would be amorphization of initially crystalline 
dust, which may occur under heavy proton bom- 
bardment due to the stellar wind. This last mecha- 
nism also explains the lack of crystalline silicates i n 



the ISM ( Carrez et al.ll2002l : iBrucato et al.ll2004 ) 



because these features resemble those observed in 
red supergiant stars, the presence of crystalline 



5.2. Features of ionized nebulae 

Ratios of fine structure lines can be compared 
with typical values found toward H II regions, 
photo-dissociation regions (PDRs), SNR, and the 
ISM. First, we have considered the ratio between 
[Si II] 34.8 /um and [S III] 33.5 /xm. Theseq two lines 
are found in the ISM, and their ratio is indicative 
of the degree of excitat ion of the gas. According 
to lSimpson et al.l (|2007l ). this line ratio is ~ 2.5 for 
the diffuse interstellar gas, and just in the range 
0.1-0.6 toward H II regions. Since the values es- 
timated for the G79. 29-1-0. 46 region are similar to 
those of H II regions (0.3 for the IRDC, 0.2 for the 
shell, and 0.5 for the central object), this suggests 
that the three components studied in this work 
have important amount of ionized gas. While this 
was expected for the shell and the central object, 
where free-free thermal emission was already de- 
tected (jffiggs et al.lll994l ). we note that the IRDC 
component is most likely picking some emission 
from the H II region DR 15, which is located only 
4' to the (south) east of the LBV star and is cov- 
ered by the IRS slit. 

Similarly, a correlation was found for the ra- 
tio between [Ne III] 15.5 /zm and [Ne II] 12.8 /im, 
and the r atio between [S IV] 10.5 Atm an d [S III] 
18.7/xm (jMartin- Hernandez etal] |2002[ ). Such 
a correlation can be used to measure the de- 
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gree of ionization of a region. For the case of 
the central object, the first ratio is around 0.7, 
and the second is around 0.8 (assuming no con- 
tamination from He II and H I at the [S IV] 
line). By placing this value in the diagram 
([Ne III I/fNe III versus fS IVl/fS Ill l) of Fig- 



ure 1 of Martin-Hernandez et al. ( 2002[ ). we have 



found that the central object lies in an interme- 
diate case between the H II regions with low and 
high ionization, and in particular it is close to 
the value of the Orion Nebula. In addition, we 
have measured the ratio between [Ne III] 15.5 ^m 
and [Ne II] 12.8 ^m for the shell and have found 
a value a round 0.3. It is interest ing to note the 
finding of IChavamian et al. I (l2Q09h that slow and 
dense shocks produce a ratio around 0.4, which 
is very similar to the value found in the shell of 
G79. 29-1-0. 46. This suggests the presence of slow 
shocks in the shell, which is completely consistent 
with the recent discovery of a slow CO shock in 
the southwestern side of the shell of G79. 29-1-0. 46 
(Paper I). 

Finally, from the ratio between [S III] 33. 4 /im 
and [S III] 18.7 /im, one can make a rough esti- 
mate of the electron density. The raticHis 1.1, 4.0, 
and 13.0 for the central object, the shell, and the 
IRDC, respectively. The ratio of ^ 1 corresponds 
to an electron den sity of ~ 500 cm" -3 (|Rubinlll989t 
Rubin et al.l Il994l ) . These values are consistent 
wit h the electron densi ty of ~ 1000 cm~^ derived 



bv lUmana et al.l ( 20091 ) for the central object of 
HRCar. For the shell and the IRDC, this ratio 
is too high to be able to derive a reliable elec- 
tron density. This could be indicative of strong 
extinction for the shell and the IRDC cases, as 
found for some posit ions in the galactic center by 
Simpson et al] (|2007l) . 



5.3. PDR signatures 

PAH emission features are present in the 

SL spectrum of the shell. This dust material 

was also detected in the LBV stars AG Car, 



^To compute this ratio, we have corrected the Une inten- 
sity at 18.7 /jm by a factor of 3. This factor for aperture 
losses apphes only to the central object, for which the flux 
densities have been calculated in the high-resolution spec- 
tra, and the line at 18.7 fim falls in a different module than 
the line at 33.4/im. For the shell and the IRDC, there is 
no need of correction because both lines fall in the same 
module. 



HR 168625, R71 ( 


Voors et al.ll2000bl: Skinneilll997; 


Voors et al. 1999 


). PAH emission is supposed to 



be formed where shocks act upon cold molec- 
ular gas. This is presumably the situation of 
G79. 29-1-0. 46, where the ionized gas is interact- 
ing with the molecular material traced by CO. 
Note that the presence of shocks is supported by 
the morphology (see Figures [U El and [71) and by 
the detection of the CO J = 4 3 line, which has 
a critical density of at least several lO'^cm"'^. 

On the other hand, the spectra of the central 
object show [Fc II] lines, whereas the same lines 
are not detected in the IRDC spectra. These two 
observational findings would indicate that the Fe 
abundance close to the star is higher than in the 
surrounding ISM. In the ISM, Fe is highly depleted 
from the gas phase because of condensation onto 
dust grains, but shocks would reprocess these dust 
grains and release Fe atoms into the gas phase. 
Therefore, a high abundance of Fe in comparison 
with the ISM would be another evidence of the 
presence of shocks close to the LBV star. 

The presence of PAHs in the spectrum of the 
shell, the detection of low-excitation fine struc- 
ture lines such as 26.0 /im [Fe II] in the central 
object spectrum, and the CO gas emission, point 
out to the existence of a PDR surrounding the 
ionized part of the nebula. Laurent et al. ( 2000l ) 



showed a typical spectrum of a star-forming re- 
gion, an H II region, and a PDR (see their Fig- 
ures 1 and 2). In our case, the spectrum of the 
shell shown in Figure |4] is an intermediate case 
between a PDR and an H II region. In addition, 
the non-detection of [S I] at 25 fim and the detec- 
tion of [S III] are cha racteristic of photon-ionized 
gas (|Micelotta et al.l[2009.) . I n fact, [S III] and 



[Si II] are typical P DR tracers (|Tappe et al.H2006t 



Umana et aLll2009l ). alo ng with [Ar III], [Ne II], 



[Ne I II ], and [S IV] (e.g. 



2002 : ICompiegne et al 



Martin-Hernandez et al 



200S 



Ghavamian et al 



20091 ). all present in the central object spectra. 

In PDRs, PAHs are non-thermally excited 
by s ingle ultraviolet (UV) photons (jSellgren 
19841 ) ■ and far-UV photons control both the ther- 
mal and chemical structur e of the neutral gas 
(jTielens fc Hollenbach|[l985l) . The high values of 
both the stellar effective temperature (Toff) and 
luminosity (L, ) favor the presence of a PDR where 
the significant UV field encounters cold molecular 
gas. However, as the CO shells are at large stellar 
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distances, the presumably high UV field may be 
diluted. 

Thus, we have made an attempt to estimate the 
incident UV field (Go) in the 24 /im inner shell. 
Go is defined as the fraction of an incident UV 
field in the ionizing UV range (i.e., between 6 and 
13.6 eV). This is usually expressed in units of the 
Habing field (1.2 x 10"'* ergcm"^ s~^ sr~^), which 
is a representative value of the mean interstellar 
UV field. Go depends on the stellar parameters as 



Go oc /uv(T'cff ) • i* • d 



(1) 



where d is the distance of the affected volume to 
the star and /uv is the fraction of the Planck func- 
tion emitted between 6 and 13.6 eV with respect 
to the total: 



/uv(reff) 



(2) 



By varying T^e between 15,000 and 30,000 K, 
and L* between 1 and 2 x 10^ L0, Go ranges from 
0.8 to 3.0 X 10^, in units of Habing field. In the 
assumed range of Teg, /uv varies from 31% to 56% 
of the total energy emitted by the star. A value 
for d of 1 pc was assumed, which corresponds to 
the distance of the 24 /im shell to the star, at a 
distance of 1.7 kpc. 

These high values of Go (of the order of 10^), 
together with moderately high H2 densities (at 
least several 10^ cm~^ due to the presence of CO 
J = 4— ^3), are similar to the conditions found in 
extreme PDRs, such as the well-know n Mon R2 
(|Rizzo et al.ll2003al: iBerne et al.ll2009l). It is also 



compa tible with the values found bv lUmana et al 



in HRCar. However, PAHs are expected 
to be not widespread but located in small regions 
where the restrictive conditions of PDRs apply. A 
zone close to the inner 2A^m shell and the CO 
J = 4 — > 3 shells is a good candidate to search for 
these fragile large molecules. 

5.4. Molecular hydrogen 

The 0-0 S(0) line of H2 has been detected in 
the three targets, while the S(l) line has only 
been identified in the central object and the shell. 
Unfortunately, we have not detected more excited 
pure rotational lines (S(2) and above), which pre- 
vent us a thorough excitation study. Despite that. 



the mere presence of molecular hydrogen in the 
close environs of a LBV star is interesting and de- 
serves further analysis. To the best of our knowl- 
edge, the only repo rted detection of H2 cl ose to 
a LBV was done bv lMorris. P.. et all (l2008l) . It is 
important to note that lUmana et al. 1 (120091) . using 
similar data sets than ours, did not found any H2 
lines in the dusty nebula HRCar. 

Vibrationally excited H2 was also discov- 
ered surrou nding part of the W olf-Rayet nebula 
NGC2359 (ISt-Louis et all Il998l ). in correlation 
with significant amounts of (probably) shocked 
CO (|Rizzo et al.ll200lU2003b^ . 

Our values of the S(1)/S(0 ) line ratio are not 
too different to those found bv lTappe et al. I (l2006h 
in a molecular cloud to an SNR, affected by both 
UV radiation and shocks. Unfortunately, a deeper 
analysis from this line ratio is not appropriate due 
to the fact that both lines correspond to different 
species, and we cannot constrain the ortho-to-para 
ratio. A better study including more lines and 
a precise spatial distribution may help to under- 
stand the nature of the emitting H2 lines, as well 
as a comparison to several PDR models (see, e . g . 
Draine fc Bertoldilll996HCompiegne et al.ll2008[) 



6. Conclusions 

In order to shed light on the distribution of 
warm dust, PAHs and CO gas in G79. 29-^0.46, we 
have analyzed infrared Spitzer data and observed 
the CO J = 2 ^- 1 and 4^-3 lines using the IRAM 
30m radio telescope and the SMT. 

An analysis of the high-sensitivity MIPS 24 fim 
image has led us to the discovery of a second outer 
dusty shell with a spherical shape, previously over- 
looked. In addition, the analysis of the SED of the 
central object has revealed the possible existence 
of a third extremely young dust shell closer to the 
LBV star, unresolved in the Spitzer images. Fur- 
thermore, the new CO data shown in this paper 
have revealed two separated gaseous shells located 
between the two dusty shells. The dynamic time 
of these gaseous components would indicate that 
the various shells have been formed during the cur- 
rent evolutionary stage. The presence of multiple 
shells around a LBV star is of great interest be- 
cause its implications for stellar evolution theories. 
This would indicate that the mass-loss history has 
been more complex than previously assumed, with 
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periods of enhanced and dimmed stellar winds in- 
stead of a continuous steady one. 

Using the IRS low-resolution spectra, we 
have compared three regions in the field of 
G79. 29-1-0. 46: the central object, the shell, and the 
IRDC. The three studied regions present a clear 
similarity at the longest wavelengths (>25/xm), 
due to the presence of a large amount of fore- 
ground/background dust at '--^ 65 K. A rich variety 
of solid-state features has been detected. Amor- 
phous silicate features in absorption are present 
in the spectra of the central object and the IRDC, 
probably due to the foreground ISM and/or to 
the presence of amorphous silicate dust in a likely 
inner nebula very close to the central LBV star. 
PAHs features in emission are present in the spec- 
trum of the shell, and a tentative detection of 
amorphous silicate in absorption is reported as 
well. Furthermore, a CO2 ice feature in absorp- 
tion is present in the spectrum of the IRDC. 

The IRS high-resolution spectra of the central 
object are full of fine structure lines, similar to 
those detected in other LBV stars. We have also 
detected pure rotational lines of II2 arising from 
optically thin quadrupole transitions in two re- 
gions of G79.29-f 0.46, the central object and the 
shell, which is the first detection of H2 in the spec- 
trum of a LBV star. 

Using different flux line ratios, we have ana- 
lyzed the physical conditions of G79. 29-1-0. 46 and 
its surroundings. In the case of the shell, indica- 
tions of slow shocks have been found. 

We have estimated an incident UV field in the 
inner dusty shell of the order of 10'* in units of 
Habing field. This high value of Go, together with 
other evidence such as moderately high H2 den- 
sities, the presence of PAHs in the shell, the de- 
tection of low-excitation fine structure lines, and 
the CO mid-J emission, suggests the existence of 
a PDR. 

This study of the multiple phase material di- 
rectly associated with G79.29+0.46 may help to 
understand its recent past and consequently the 
evolutionary mechanisms of this kind of object. 
In this sense, G79.29-f0.46 becomes a paradigm 
of LBV nebula and deserves further studies. In 
particular, a follow-up study of complex molecules 
should also be carried out, because they can pro- 
vide information about the physical parameters 



in detail. In addition, similar studies to the one 
shown in this paper for other LBVs and evolved 
massive stars will also help to understand this 
short, but crucial, evolutionary phase of high-mass 
stars. 
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